Divergence of the entanglement range in low dimensional quantum systems 
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We study the pairwise entanglement close to separable ground states of a class of one dimensional 
quantum spin models. At T=0 we find that such ground states separate regions, in the space of the 
Hamiltonian parameters, which are characterized by qualitatively different types of entanglement, 
namely parallel and antiparallel entanglement; we further demonstrate that the range of the Con- 
currence diverges while approaching separable ground states, therefore evidencing that such states, 
with uncorrelated fluctuations, are reached by a long range reshuffling of the entanglement. We gen- 
eralize our results to the analysis of quantum phase transitions occurring in bosonic and fermionic 
systems. Finally, the effects of finite temperature are considered: At T>0 we evidence the existence 
of a region where no pairwise entanglement survives, so that entanglement, if present, is genuinely 
multipartite. 

PACS numbers: 03.67.Mn, 75.10.Jm, 73.43.Nq, 05.30.-d 



I. INTRODUCTION 

Quantum fluctuations may disorder the ground state 
of a system, especially at low dimensions. Paradigmatic 
example in this sense are quantum phase transitions [l|, 
where different phases can be achieved at T — by ad- 
justing a control parameter of the system. Paradoxically 
enough, quantum effects can provide also classical-like 
ground states (CGS). In fact, for certain values of the 
control-parameter, quantum fluctuations may become 
completely uncorrelated; thereby the ground state of the 
system gets factorized and identical to the lowest-energy 
state of the classical counterpart of the original quantum 
system. This phenomenon was evidenced by Kurmann 
et al. in the early eighties, for 5=1/2 Heisenberg an- 
tiferromagnetic chains in an external magnetic field h. 

The study of entanglement in quantum many-body 
systems has beenproviding a new angle in statistical 
mechanics 0, 0, H IE 0> particularly at low temper- 
ature where cooperative phenomena are dominated by 
quantum mechanics. Thanks to a proper analysis of cer- 
tain entanglement estimators, the result by Kurmann et 
al. was recently retrieved _|8| and generalized to two- 
dimensional spin systems moreover, numerical evi- 
dences arose for it to hold in spin chains with long-ranged 
interaction [To| . 

Special interest has been devoted to bipartite entangle- 
ment of formation in connection with quantum critical- 
ity: In fact, though quantum critical points are rather 
marked by the enhancement of multipartite entangle- 
ment [HllEl> the variation of the pairwise entanglement 
at criticality captures the non analyticity of the ground 
state of the system [l2j]- On the other hand, the naive 
guess that the range of pairwise entanglement should di- 
verge at a quantum phase transition, in analogy with the 



divergence of the correlation length of the two-point cor- 
relators, has never been evidenced [1,0. 

In this paper we show that, in the space of the Hamil- 
tonian parameters, the special points where CGS occur 
(hereafter called separable points) mark a separation be- 
tween regions characterized by different types of entan- 
glement, (called antiparallel and parallel entanglement 
in Ref. [13), which correspond to qualitatively different 
spin configurations. The transition from one region to 
the other is found to be characterized by the divergence 
of the range of pairwise entanglement in the immediate 
neighborhoods of the CGS (see Eqs. and (|T3jl below). 

Spin-off in systems of strongly interacting bosons are 
also found: We evidence that the superfluid-insulator 
quantum phase transition at commensurate filling is in 
fact a transition between a phase (superfluid) where 
solely particle-hole entanglement is present and a phase 
(insulator) with no pairwise entanglement at all: the 
particle-hole entanglement is easily seen to correspond 
to antiparallel entanglement and the range of the concur- 
rence is found to diverge while approaching the transition 
(see the paragraph below Eq. lj*P3|) ). 

Finally, we study how robust CGS are with respect 
to temperature: We evidence the emergence of a region 
in the h — T plane, fanning out from separable points 
(see Fig. |3J) where pairwise entanglement vanishes. In 
such region, if entanglement is present in the system, 
it necessarily is of multipartite type. Our study does 
also shed light on the result by Kurmann et al. (see the 
concluding paragraph). 
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II. MODELS 

We focus our attention on the class of one-dimensional 
spin models described by the Hamiltonian 

^(jxiii/jiz) = J (jx sf sf + 1 +j y sf S- + 1 +j z Sf s* +1 

i 

- hs;) , (i) 

where i runs over the sites of the chain, Sf [a = x, y, z) 
are quantum angular momentum operators correspond- 
ing to S — 1/2, j a are the anisotropy parameters (\j a \ < 
1), h = gfi^H/J is the reduced magnetic field, and J is 
the exchange integral, assumed positive and hereafter set 
to unity. 

In Ref0 it was demonstrated that CGS are obtained 
for h = hf = yj{j x +j z )(j y +j z ); although the model 
(JTJ cannot be solved exactly for generic j a , the above 
result is rigorous. In order to analyze quantum corre- 
lations, which are crucial for understanding the behav- 
ior of the system at and near a separable point, we re- 
strict the parameters in the Hamiltonian (JTJ so as to 
rely on exact results: We therefore consider the solvable 
cases W(l + 7, 1—7,0) = Hxy with < 7 < 1, and 
TC(1, l,jz) = Ti-xxz corresponding to the XY and XX Z 
models in a transverse field, respectively. 

The quantitative analysis of the pairwise entanglement 
between two spins sitting on sites I and m of the chain 
is addressed via the concurrence C r with r = \l — m\ 
(translational invariance is assumed) 14]. In terms of 
correlation functions g" a = (SfSf.J) and magnetization 
M z — (S z ), the concurrence reads [l5| 

C r = 2max{0,C;,C;'} , (2) 

C' r = + , (3) 

C' r ' = \9 X r X -9 V r V \+9 Z r Z -\, (4) 

where C' r and C" measure the pairwise entanglement 
related with the occurrence of antiparallel and parallel 
Bell states, respectively, as discussed in Refill both for 
pure and mixed states. We will also consider the one- 
tangle Em 

r 1 = l-AY J Ml, (5) 

a 

the two-tangle 

r 2 = 2^C r 2 , (6) 

r 

and the ratio Tijr\ which estimates the fraction of the 
total entanglement stored in pairwise correlations [S| . Al- 
though gf a and M z do not show any anomaly at a sepa- 
rable point, they unveil it when combined in C r , which is 
found to drop to zero in a non-analytic way at this point. 
Such circumstance comes together with the vanishing of 
Ti, and with a highly non-trivial behavior of the ratio 
r 2 /ri i. 



III. RESULTS 

A. Zero Temperature 

Let us first consider the completely integrable 0, El 
XY model in a transverse field. Besides the quantum 
critical point h = h c = 1, its T = phase diagram is 
characterized by the circle h 2 + 7 2 = 1 (hereafter called 
the separable circle) where CGS occur in the form 

\GS XY > = U\4>? Y > , 

i 

\tf Y > = (-iy CO s( d -i)\u>+sm( e f)\u> , 
cos( ^ } = VIT^ ' 

where \<fif Y > is the state of the spin sitting on the z-th 
site. 

Since the early papers on the model it is known that 
the functional form of M z and gf a depends substantially 
on whether the parameters of the system fix a point sit- 
ting inside, outside, or on the circle itself [ljj; however, 
it had never been noticed that the simplicity of the ex- 
act solution at hf Y = y/l — j 2 is ultimately due to the 
factorization of the ground state. In fact, this is clearly 
evidenced by the concurrence, whose expression on the 
separable circle reads 

C r = 2C' r = 2C'; = j-t- + 2M 2 - 1 - Vr , (7) 

and hence, being M z = — 7)/(l + 7) , 

C r = C'r = C'r' = , Vr . (8) 

Moreover, it is C' r > (and C" < 0) inside the circle, 
and C' r < (and C" > 0) outside the circle, no matter 
the sign of the exchange interaction and the value of r. 
According to the analysis presented in Rcf 13 this means 
that inside (outside) the circle pairwise entanglement ex- 
clusively originates from the occurrence of antiparallel 
(parallel) Bell states. 

We remark that whether the system has parallel or 
antiparallel pairwise entanglement at T = uniquely de- 
pends on the Hamiltonian parameters 7 and h: As a 
consequence, we can draw an "entanglement phase di- 
agram" in the h — 7 plane, where different phases are 
characterized by the presence of parallel or antiparallel 
entanglement. The separable circle h 2 + 7 2 — 1 marks 
a boundary in such diagram (see Fig. ^) suggesting that 
the occurrence of a CGS is a necessary step for switch- 
ing from parallel to antiparallel entanglement. We notice 
that the same scenario emerges from the numerical anal- 
ysis of more complicated models, both in one and two 
dimensions 

As the transition from parallel to antiparallel entan- 
glement involves the whole system, we study how the 
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FIG. 1: Entanglement phase diagram in the h — 7 plane. 
Points on the circle correspond to models whose ground state 
is separable. 

pairwise entanglement propagates along the chain in the 
vicinity of a separable point: for doing that, we introduce 
the range R of the T — pairwise entanglement, which 
is defined as the maximum distance between two spins 
such that the concurrence is non vanishing: 

C r > for r < R and C r = for r > R . (9) 

We underline that the exact vanishing of C r for r > R 
and h h{ follows from the definition of the concurrence 
Eq. @, in the sense that C r vanishes whenever C' r and 
C" are both negative. On the other hand, at the factor- 
izing field, C r = C' r = C'l = for all values of r due 
to the fact that the correlation functions do not depend 
upon r on the separable circle. In Fig. [3 we show the 
exact results for C r with r ranging from 1 to 5. Results 
for larger r are also available and show the same quali- 
tative behaviour: C r fans out from the separable point 
with non-zero derivative, reaches a maximum, and then 
vanishes, both for h > hf and h < hf, though not sim- 
metrically with respect to hf. 

The observed behaviour suggests the divergence of R 
for h — > hf\ by expanding the exact expressions of the 
correlation functions up the first order in (h — hf) we 
find [13 

Cr = ^Ll\h - h { \ + o((h ^ hf) 2 ) , (10) 



where T = - j)/(l + 7). Eq. (UUJ) confirms that 
all C r get progressively positive for h approaching the 
factorizing field, as clearly seen in Fig. [21 This means 
that the range of the concurrence R diverges at hf. 

In order to analyze how R diverges with the field, we 
push forward the expansion in (h — hf), meanwhile con- 
sidering the large- R expressions for the correlation func- 
tions, given in Refs^J and |^. The result for h > hf 
reads 

C'l = ^Ll{ h -hf)-[A 2 + B{r)]{h~hff+0{{h-hff) 
47 

(11) 



FIG. 2: Concurrence for the XY model at T = 0: curves are 
C r versus h for r = 1, 2, 3, 4, 5 (top to bottom), and 7 = 0.5. 
The inset is a zoom near the factorizing field, in logarithmic 
scale. 

where A 2 = T 2 (3 + 7)/(327 3 ) and B(r) is a coefficient 
which vanishes for r — > 00. The range of the con- 
currence is implicitely obtained from Eq. (jl 1|> by re- 
quiring = 0. In fact, for sufficiently large r, it is 
(A 2 + B(r))- 1 ~ A- 2 (l - B(r)/A 2 ), and hence, by sub- 
stituting into Eq. (|ll(l . C" is found to vanish both for 
h = hf and h - h { — r 2r_1 /(47/I 2 ), leading to the loga- 
rithmic divergence 

R XY cx ^ln 1^1 \ \n\h-hf\~\ (12) 

where we have introduced the symbol R XY to make clear 
that the functional form of the divergence is in general 
model-dependent. 

For h < hf the expression for C' r is different from 
Eq. Qllfl. and to this difference we ascribe the asymmet- 
ric behaviour of C r with respect to the separable point 
observed in Fig. [2 however, for h — > hj the above result 
is still valid, though just for the (thermal) ground state 
with unbroken symmetry 0,H2- It is to be noticed that, 
in the thermodynamic limit (here considered) and while 
approaching the CGS, the fact that the concurrences Cf 
become progressively positive for larger and larger r goes 
together with their getting vanishingly small: this is due 
both to the monogamy of the entanglement [23L 0] and 
to the proximity of the separable point itself. 

The divergence of R implies, as a consequence of the 
monogamy of the entanglement, that the role of pairwise 
entanglement is enhanced while approaching the separa- 
ble point; in fact, the ratio t%It\ is found to have a cusp 
minimum at the critical point and to increase while mov- 
ing towards the CGS, in full analogy with the result of 
Refs@ and In the particular case of the Ising model 
(i.e. 7 = 1), we find that for h — * hf = the ratio 
Tijr\ goes to unity at the separable point. For 7 ^ 1 
and hf < h < h c our data show that t%It\ monotoni- 
cally increases for h — > hf and that the value (t2/ti)L+ 
increases with 7. 
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We remark that the divergence of R XY while approach- 
ing the separable circle cannot be recognized as a critical 
effect, since the ground state is non-singular at hf and 
the long-ranged pairwise entanglement does not survive 
neither inside nor outside the separable circle. 

A complementary view on the physics of CGS is ob- 
tained by the analysis of Hxxz, that can be done re- 
sorting to Bethe ansatz results |25|. In this case hf 
coincides with the saturation field h s = (1 + j z ), and 
\GS XXZ >= I1J Tt>- Also > and distinctively from the 
XY case, the factorized state extends over a finite por- 
tion of the h — j z phase diagram of the model. In fact, 
h s separates a gapless quasi-ordered phase (with power 
law decay of the in-plane correlation functions) from the 
gapped, fully polarized phase (with (S z ) = 1/2). Due to 
the in-plane symmetry of the model (implying gf x = g^ v ) 
C" is always negative and therefore pairwise entangle- 
ment, if present, is of antiparallel type. Combining the 
exact results of Refs|2^ andl27l we find that R diverges 
approaching the band transition as 



R 



xxz 



cx (h s — h) 



(13) 



where 



1 



7rtan(^?) tan(7rry) 



and r\ = — arcos(— j z ). 



(14) 

The divergence of R in the isotropic case (j z — 1), specif- 
ically determined in Refll3l results from Eq. (|13fl with 
8 = 2. We underline that saturation is not related to 
a spontaneous symmetry breaking, and the divergence 
of R while approaching h s does not mark a critical phe- 
nomenon. 

Let us now consider strongly interacting hard-core- 
bosons/spinless-fermions in Id: We shall find that the 
occurence of CGS plays a fundamental role for such sys- 
tems. 

Hard-core bosons with repulsive Coulomb interaction 
are described by the quantum lattice gas |28|, whose 
dynamics is described by 'Hxxz phrasing the spins in 
terms of hard-core bosonic operators: a = S~ , a 1 = S + , 
a^a = S z + 1/2. The relevant energies t-^j x (here 
jx=jy=l), U^j Zl and + j z are the hopping am- 

plitude, the Coulomb interaction, and the chemical po- 
tential, respectively. By this mapping, the superfluid and 
insulating behaviors of the quantum lattice gas at com- 
mensurate filling are related to the quasi-ordered (par- 
tially filled band) and fully-polarized (filled band) phase 
of the XXZ spin model, respectively; therefore, the band 
transition observed at [i=t+2U, corresponds to satura- 
tion occurring at h s =(l+j z ). 

Based on the above analysis, we state that the 
insulator-superfluid band transition is characterized by 
the divergence of the range of the concurrence. Remark- 
ably, the antiparallel character of the pairwise entangle- 
ment present in the XXZ model reflects the fact that 
close to the superfluid-insulator transition exclusively 
particle-hole entanglement plays a role. Arguments along 




FIG. 3: Contour plot of ti in the h — T plane, for 7 = 0.3 
(i.e. hf = 0.9539...). The white area indicates where T2 = 0. 
Condition l|15|l is fullfilled below the dashed line, which is 
defined by (H) = E scp . 



the same line can be applied to spinless fermions mod- 
els obtained via Jordan- Wigner transformation of Hxxz 
[23 : The band-transition there observed is from an insu- 
lating regime to a gapless phase equivalent to a Luttinger 
liquid. 



B. Finite Temperature 

We now switch on a finite temperature in the system. 
We consider questions like: What is the effect fanning 
out from CGS on the thermal (mixed) states of the sys- 
tem? Particularly: How meaningful is the characteri- 
zation of the system in terms of parallel or antiparallel 
pairwise entanglement for mixed states? To answer these 
questions we consider the XY model where both paral- 
lel and antiparallel entanglement are present at T = 0. 
The analysis of t-i evidences that in the h — T plane it 
exists a region, fanning out from the CGS, where pair- 
wise entanglement vanishes (white region in Fig. [2} and 
the entanglement, if present, is shared between three or 
more parties. This means that a CGS may evolve into 
a quantum mixed state with genuinely multipartite en- 
tanglement by increasing temperature. In order to deter- 
mine whether this happens or not, we need to know if 
entanglement is present in the system: the one-tangle, 
which accomplishes this task at T = 0, is not a proper 
estimator for the entanglement content of the system at 
finite temperature; therefore, we have to refer to some en- 
tanglement witness. Following the results by Toth [30| . 
entanglement is present if 



(H) 



E scp < 



(15) 



where E scp is the ground state energy of the corrispond- 
ing classical model. The region below the dashed line in 
Fig-EDis where condition ll5l is fulfilled, i.e. where entan- 
glement of whatever type is present in the system. We 
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FIG. 4: C r versus r for 7 = 0.01 and ft = 1.1. Circles, 
squares and triangles correspond to three different tempera- 
tures: Ti=410" 3 , T 2 =4.710~ 3 , and T 3 =510" 3 , respectively. 
Full symbols mean C r ~C r and empty symbols mean C r =C". 



further observe that, in contrast to the analysis of the 
ground state, at finite temperature we cannot character- 
ize the two separate phases of parallel and antiparallel 
entanglement. In fact, the two types of entanglement 
(though well defined also for mixed states) can swap 
by varying T and/or r (see Fig. QJ. The exchange be- 
tween parallel and antiparallel entanglement occurs in a 
non trivial way, that ultimately produces the tempera- 
ture "reentrance" of seen in Fig [31 We also find a 
regime where C r can be a non monotonic function of 
r, so that, for instance, one spin is not entangled with 
its nearest neighbor while being entangled with its next- 
nearest neighbor. Such situations occur in the vicinity of 
the region where C r vanishes, as seen in Fig. 21 and it is 
due to the non-monotonic behavior|l9j of the correlation 
functions. 



IV. CONCLUSIONS AND PERSPECTIVES 

Summarizing, we have studied the occurrence of CGS 
in relation with pairwise entanglement, in a class of 
5=1/2 spin chains. Our results show that at T = 
the space of the Hamiltonian parameters is divided into 
regions where either exclusively parallel or exclusively an- 
tiparallel pairwise entanglement is present, no matter the 
distance between the considered spins: Therefore, an en- 
tanglement phase diagram can be unambiguously drawn. 
Transition lines in such diagram corresponds to separable 



ground states, and are further characterized by the fact 
that the range of the concurrence diverges while moving 
towards them. Due to the monogamy of the entangle- 
ment, such divergence goes together with the asymptotic 
vanishing of the value of the concurrence itself. 

We further provide (to our knowledge for the first time) 
an explanation of the phenomenon described by Kur- 
mann et aZ. 2?] : the factorization of the ground state is 
a necessary step for antiparallel entanglement to be fully 
replaced by parallel entanglement. We observe that in 
the global reshuffling of the ground state which leads to a 
CGS, and involves all the spins of the chain, a long range 
entanglement appears as a crucial ingredient. Moreover, 
our results for the entanglement ratio ti jr\ confirm that 
multipartite entanglement dominates at a genuine quan- 
tum critical point, while pairwise entanglement is essen- 
tial for understanding the mechanism leading to CGS in 
quantum systems. 

Our analysis is of relevance also for bosonic systems un- 
dergoing a supcrfluid-insulator transition. It is intriguing 
to conjecture that the divergence of the range of C r at a 
CGS goes beyond model dependency. In fact, Anfossi et 
al. [Tj have observed a similar divergence of the range 
of C r also in the bond-charge extended Hubbard model 
at certain transition lines. 

For finite temperature the entanglement in the system 
cannot be characterized by the single parameter h and 
a much more complicated scenario emerges: In particu- 
lar we find that, by increasing T, the factorized (pure) 
ground state evolves into a thermal (mixed) state with 
null pairwise entanglement: this opens the possibility 
for the existence of an experimentally accessible finite- 
temperature region where entanglement, if present, is 
genuinely multipartite. 

Finally, we notice that the possibility of controlling via 
a proper tuning of the external magnetic field whether 
two spins are entangled or not, and whether they share 
parallel or antiparallel entanglement, might be of inter- 
est both from the experimental point of view and for 
applicative purposes. 
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